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Abstract 

We define a hierarchy of dynamic relaxed gas spheres as solutions 
^- ■ of the Poisson equation coupled to a hierarchy of approximations of the 

l/~) ' Liouville equation leading, when this equation is satisfied, to the well- 

ed . known isothermal gas spheres in the static case, but also to a new form 

of dynamic maximum relaxation. The two previous steps of the hierarchy 
correspond to an increasing degree of local relaxation at the center of the 
configuration. 
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1 Introduction 

i—t , Embodied into the concept of equilibrium there are two ingredients. The first 

P*J' and most primitive is the idea of stillness. The second one is the idea that 

equilibrium is the latest stage of a macroscopic process to which a system has 
been driven by a microscopic mechanism at work in every cell of the system. It 
would be superfluous to distinguish these two ingredients if we assumed that the 
process under consideration necessarily leads to stillness, but we do not need to 
do that. Here we distinguish between the concept of equilibrium proper with 
the usual connotation of stillness attached to it, and the concept of maximum 
relaxation which is a concept that describes some terminal dynamic stages of 
evolution. 

We have already considered this distinction beforejj in the framework of 
General relativity and more precisely in the description of the matter contend 
of some Robertson- Walker models of universe using Relativistic kinetic theory. 
Since the Universe is not stationary, i.e, there does not exist any global time- 
like Killing vector-field, but expanding it is not possible to assume that the 
cosmological fluid is in equilibrium in the usual sense. But it is possible to 
imagine that this fluid is in an state of maximum relaxation, which in the 
quoted reference we called generalized equilibrium. We believe now that the 
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concept of maximum relaxation may be useful also in more classical contexts, 
either to describe laboratory systems with time-dependent boundary conditions 
or external perturbations, or in astrophysical descriptions of self-gravitating 
systems in the strict framework of Newtonian gravity. This is more specifically 
the problem that we discuss in this paper. We hope that this restriction to a 
more classical problem may help to uncover the usefulness of the concept of 
maximum relaxation. 

Actually the description in a particular case of an state of maximum re- 
laxation is not the unique goal that we pursue in this paper. We are equally 
interested in describing some intermediate steps that may settled down in sys- 
tems which do not reach uniformly maximum relaxation. It might well be that 
some of these partial stages of relaxation become interesting by themselves. 

It is possible to distinguish three domains in the physics of the kinetic theory 
of gases: i) the physics of equilibrium, ii) the physics near equilibrium where 
perturbation theory can be used, and iii) the jungle where the most difficult 
problems arise. This paper is a modest contribution to a partial clarification of 
a small territory of this non perturbative jungle. 

2 Particular field integrals 

Let us consider a spherical mass distribution with Newtonian gravitational po- 
tential V{t, r) and a massive particle free-falling in the corresponding gravita- 
tional field according to the Newtonian equations of motion: 

— = -^d J V{t,r), v* = —, r=^7^r, e t = l (1) 

We shall say that the following quadratic expression: 

1 
2 

is a spherically symmetric Field integral if: 

• Aij and Bi have the following tensor structure: 

Aij = M6ij + Nmrij, Bi = Pm (3) 

with: 

Hi = Sijn 3 , ri = r~ x 1 , (4) 

M, N, P and C being functions of t and r. 
The total derivative of I with respect to t is independent of the velocity: 

dl „, 

Tt =F(t,r) (5) 

where F is a function that we shall call the source of the field integral. 



i = -A t y v j + By + c (2) 
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Requiring / to satisfy this last equation taking into account 0, yields, to 
start with, the following set of conditions: 

8 {k A lf) = (6) 

^dtAij + d {l B 3) = (7) 

- d^dkVAij + d t B % + d t C = (8) 

-6 jk d 3 VB k + d t C = F (9) 

And taking into account |3J we obtain the following system of equations: 

M' + =0, N' =0 (10) 

r r 

2P ■ 2P 

M+ — = 0, N + 2P' =0 (11) 

r r 

-V'{M + N)+P + C =0 (12) 

- V'P + C = F (13) 

where a dot means partial derivative with respect to t and a prime a partial 
derivative with respect to r. 

Integrating the system of Eqs. EQ-U3 we get: 



N = -vr\ M = vr 2 +n (14) 

P = -\iir (15) 

C = fiV + -hr 2 + a (16) 

where v is a constant, and fi(t) and a(t) are two functions of t. Taking into 
account [15] and [l(j Eq. yj] becomes: 

where /i' 3 ' means the third derivative of fi with respect to t. 

The concept of a field integral is an obvious generalization of the concept of 
constant of motion to which it reduces if F = 0. This concept can be useful 
only if some other considerations restrict its source F. This will be the case in 
the sequel of this paper. 

Collecting the above results, I can be written as: 



I = E + ^J 2 (18) 

where: 

E = p{-v 2 + V)- -prvn+-flr 2 +a, J 2 = r 2 [v 2 - [vnf]. (19) 

J, which is the angular momentum per unit mass with respect to the center of 
the mass distribution, is a constant of the motion, but E, which is proportional 
to a generalization of the energy per unit mass to which it reduces if p = 1 and 
V = is a field integral: 

d ^=F, ^=0 (20) 

dt ' dt y ' 

3 Coupling gravitation to a kinetic equation 

We assume in this paper that the density p(t, x l ) of the gas under consideration 
is given as the mean value of a kinetic distribution function f(f,x % ,tP). More 
precisely we assume that the density p in the Poisson equation 

AV = p (21) 

where units have been chosen such that: 

4ttG = 1 (22) 

is given by: 

p(t,x i ) = f f{t 7 x i ,v J )dv 1 dv 2 dv 3 (23) 

where / > is a solution of a kinetic equation: 

d t f + v l dJ~5^d l V^- = B (24) 

with B(t,x z ) being a function independent of the velocity that has to be zero 
when equilibrium or maximum relaxation is reached. 

Since we are interested only in spherically symmetric configurations, B, p 
and V will be functions of t and r, the equation |2l] becoming: 

V" + -V = p, (25) 

r 

and / will be a function of t, r, v 2 and vri. 

Finally, we shall assume that the domain on which O and 24 apply includes 
the center of the gas sphere and is bounded by a sphere of finite radius R. And 
also that in this domain the density and the mass M are finite. Consistently 
with these assumptions we shall assume without any loss of generality that 



V(t, 0) = 0, but we have to keep in mind that the value of this quantity should be 
chosen appropriately to match the interior solution to the solution corresponding 
to a next layer beyond the surface r = R. Complete models and the matching 
problems are not considered in this paper which deals only with the central 
cores of gas spheres where presumably some degree of relaxation has already 
been achieved. 

4 The road to equilibrium or maximum dynamic 
relaxation 

We assume from now on that / is a distribution function that has the following 
form: 

f = mexp(-I) (26) 

where m is the mass of the particles of the gas and / is a field integral of the 
type that we considered in section 2, with N, M, P and C as in [lj, [l5| and [16], 
and with F being given by Eq. [17]. Therefore the r-h-s of p4| is: 

B = -fF (27) 

Evaluating the integral 03 we obtain: 

/x J /^(l + /i L vr z ) 

We follow below the hierarchy of steps leading to equilibrium or to maximum 
relaxation described by an ordered set of restrictive assumptions on the function 
F(t, r). We shall consider a power expansion of the this function with respect 
to the variable r: 



F(t, r) = F (t) + Fi(t)r + ■■■ + ^F n (t)r n + ■■■ (29) 

and we shall examine the models which are obtained assuming that more and 
more terms, starting with Fq, of this expansion are zero. 

Concomitantly with E3 we shall consider the power series expansions of 
V (t, r) and p(t, r) derived from [28f 

V(t,r) = V (t) + Vi(t)r + ■■■ + \v n {t)r n + ■■■ (30) 
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p(t,r) = p (t) + Pl {t)r + ■■■ + ^Pn(t)r n + ■■■ (31) 

Integrating \23 term by term we obtain: 



71 — 1 

Vi = 0, V n = — — p n _ 2 , n > 1 (32) 

n + 1 



the first of these equalities coming from the assumption that the density is finite 
at the center. 

First step. Assuming that Fo(t) = in Eq. p"7J yields: 

jjlVq + /jVq + & = 0, or (iVa + cr = Constanti, (33) 

and therefore using ES and p2 we can write: 

V-2 = —po = aft ' , with a = Constant2- (34) 

o 

From Vi = it follows that F\ = 0. Also from |8| and Vi = it follows that 
pi = and therefore from [32] we have also V3 = 0. 

Second step. The next step is then to assume also that F% = 0, i.e.: 

pV 2 + \pV 2 + ^ = (35) 

and using Hj we get: 

/i (3) + ap- 3/2 p = (36) 

This is the first relevant result of this paper as it gives the equation that governs 
the dynamics of the central core of the spherical gas configuration. 
This equation has two types of solutions: 

• A. Those with (i = 

• B. Those with fi non identically zero 
Type A. In the first case Eq. O becomes: 

fiV + a = F (37) 

therefore from p3 it follows that Fq = 0. In fact, using repeatedly Eqs. p8|, 
p2| and [37] it is easy to prove that F — 0, meaning that / is a solution of the 
Liouville equation, and that p — and V n = for n > 0. The potential V(r) 
will be obtained integrating the non-linear differential equation: 

V" + V = 3ap-^ C ^-f\ (38) 

r 1 + (J, ^vr l 

If v = then the configurations thus obtained are the well-known static 
isothermal gas spheres with temperature: 

where k is Boltzmann's constant. 

Type Bl From now on we assume that we are dealing with the case B above, 
i.e. p 7^ 0. In this case the potential V(t,r) is a solution of the following 
equation: 



V" + V = 3a M - 3/2 Pl ^ UM rT "' ' (40) 



r 1 + p- 1 vr 1 

Type B2 Since from ^ and p\ = it follows that V3 = and therefore also 
F3 = the next step consists in requiring F4 — 0. This yields: 



3/iy 4 + jutfi = (41) 

From B3 and EM we have: 

y 4 = ^/,- 3 / 2 (-a M - 1/2 - lii + J/i-V - 2„- V), (42) 

Deriving this expression and substituting in f4^fl we obtain: 

/i (3) + rjtT 1 /^ + (2a/!~ 3/2 - tM~V + 2i//i~ 2 )/i = (43) 

and using in this equation [36| we get at this step the following equation: 

jl = -lap- 112 + ^"V - 4/x" V (44) 

Going back to [42] and |3^ we see that this implies that: 

Vi = 0, p 2 = (45) 

From £8] and the preceding results we have p% = and therefore also V5 = 
and F 5 = 0. 

The equation to integrate to get the potential V is: 



y, + 2 y , = 3aM . 3/2 exp(-^ + (l fl M- 1 / 2 + .- 1 ^ 2 ) 
r 1 + /i i ^)r 

Type 55. This will be the last step. From |28| and p2| we have: 



(46) 



Ve = HV M ~ 7/2 (47) 

Requiring F§ — we obtain: 

4/iy 6 +nV 6 =0 (48) 

and using ^ we get: 

ap- 7/2 pv = (49) 

and since we are considering here only dynamic configurations it follows that: 

v = (50) 

The evolution equation [|| becomes in this case: 



/2 + 2a/i- 1/2 -i^V 2 = (51) 

This is the case of maximum relaxation since, as can easily be proven, the results 
already obtained imply that: 

p{t)=p Q (t) = 3a^ 2 , V(t,r) = V a {t) + ^a^ 2 (t)r 2 , F(t,r) = (52) 

Vo(t) is at this stage an arbitrary function of t to be fixed by matching the core 
solution to its exterior gravitational field. 

5 Explicit solutions 

We are in this section interested on the explicit solutions of Eqs. 36, E3 and 
pl| , corresponding to the first and second stages of dynamic partial and total 
relaxation in the domain where the variable /i is positive. 
Eq. pra has an obvious first integral. Namely: 

b = jl - 2a/i~ 1/2 (53) 

Considering t as a function of jjl an elementary process of integration leads to: 



(i = J2bn + 8a-y/ji+2ci (54) 

from where we obtain the general solution of |3q: 

t - 1 = f^ dy = (55) 

Jii(to) y/Zby + 8ay/y + 2ci 

the three constants of integration being no, b, and c\ which are related to p(to), 
/i(io) and /i(io) by the formulas: 

b = p{t )-2ap{t Q )- 1 ' 2 , Cl = i/i(t ) 2 - ju(to)/«(to) ~ 2 flv /^) (56) 

The roots of the radical in Eq. |3 are: 
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^Ti ± = -{-2a± y/Aa^-bd) (57) 

and the values of jl for these values are: 

V4a 2 - bci 
H± = ± — . (58) 

Therefore if -JJi + are real and positive then /i + corresponds to a minimum of 
the function /i and /i_ corresponds to a maximum. 

From the preceding considerations it follows that the solutions of Bq can be 
classified according to the following types: 



• Type II: ci > 0, b > 0. The function /i has neither a maximum 

nor a minimum in the domain of interest (0 < /x < oo). If the initial value 
of jjlq is positive fj, will increase without limit; otherwise it will reach the 
value zero. 

• Type 12: c± > 0, b < 0. The function (i has a maximum, /i_. 

• Type III: c\ = 0, 6 > 0. The function /x has neither a maximum nor a 
minimum. 

• Type 112: c\ = 0, 6 < 0. The function /i has a maximum, /x_. 

• Type III1: c\ < 0, 6 > 0. The function /x has a minimum, /x + . 

• Type III2: c\ < 0, 6 = 0. The function /i has a minimum, /x + . 

• Type III3: C\ < 0, 6 < 0. The function /x has a maximum, /x_, a minimum 
/i + , and is periodic; the period being: 

Period = 2 f dy (59) 

Let us consider the following function of t: 

S(t) = }}, + 2atr 1/2 - |m _1 A 2 + 4 Ai _1 ^ (60) 

so that Eg. [44| can be written as S — 0. Evaluating the derivative of n(J and 



using Eq. p& we readily get: 

S , = - / u _1 AS'. (61) 

This means that S — is a conditional first integral of Eq. f36f if a solution 
li{t) of this equation is such that S(to) = for a particular value of t then it 
is zero for any time. Equivalently, the solutions of [44| arc the solutions |36| for 
which the initial conditions have been constrained to satisfy [ll|. 
Using pa and p4 in Eq. pq yields: 

C! = 4i/ (62) 

Therefore the general solution of Eq. E3 can be described as the general solution 
of pq given above requiring C\ to satisfy p2. In particular from pG it follows 
that the solutions of type II above correspond to configurations of maximum 
relaxation. 



6 Examples 

Let us introduce the temperature function T(t) as the following obvious gener- 
alization of the temperature parameter B3: 

r w - m (63) 

and let us consider a system of units such that besides E2 one has: 

Po {t Q ) = 3, nfo) = 1 (64) 

where to is some particular value of t. The units of time (T), length (L) and 
mass (M) of this system measured in the MKS system of units are as follows: 



f=(^irGp (to))- 1/2 s, L = Tfi(toy 1/2 m, M = (AirC)' 1 L 3 T~ 2 kg (65) 

Each set of values of the parameters b, c\ and the total mass M will de- 
fine a three parameter family of physical examples depending on the assumed 
physical values of po(io)) A*(io) an d M measured in the MKS system of units. 
Equivalently T(to) can be used if we know m. Here we shall assume that m is 
the mass of an hydrogen atom. 

We are going to discuss two particular classes of examples belonging to two 
different types among those discussed in the preceding section. 

Our first example assumes that: 

b =-2.2, ci = -1.8 (66) 

This corresponds to a type III3 and therefore fi(t) will be a periodic function. 
Assuming that to is a time corresponding to a maximum value of 
fj,, then from B4 and 67] it follows that: 

o = l, /u_ = 1, ii + = 0.8 (67) 

and performing the integral M we obtain: 

Period = 5.4 (68) 

From g5| we have that the radius R and the mass M of the gas configuration at 
any time are related by the formula: 

M = AirR 2 {t)V'(t,R(t)) (69) 

the potential function V(t, r) being obtained by integration of Eq. Eo[ Assum- 
ing, for example, that M — 7.2 and v = a numerical integration yields the 
following maximum and minimum values of R: 
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The maximum value of the central density corresponding to the assumptions 
made above is: 

{po)max = Po(to + Period/2) = 5.5 (71) 

Remark.- In this example the total mass M compatible with the values |66| 
has to be less than w 13. 

To see a physical model corresponding to the above parameters we have to 
choose a value of the central density at time to and also a value at the same 
time for the temperature. Let us assume as an example that: 

p (t ) = 0.5kg/m 3 Temp(t ) = 10 7 K (72) 

This corresponds to the following values: 

T = 8.5xl0 4 s, Z = 4.2xl0 8 m, M = 1.2xl0 25 kg (73) 

From [38|, |7(J and [7l] we have then: 

Period = 5.4 x T s = 4.6 x 10 5 s (po)max = 0.9 kg/m 3 (74) 

and: 

R min = 2.5 x L m = 1.1 x 10 9 m, R max = 2.8 x L m = 1.2 x 10 9 m, (75) 

Mass = 7.2 x M = 9.0 x 10 25 kg (76) 

Our second example assumes that: 

b = -4, ci = (77) 

This corresponds to a case of maximum relaxation of type 112, and most of its 
physical behavior is explicit in|52j and |55|. In this case also the function fi(t) has 
a maximum (T(t) has minimum). Let us assume that to is the value of the time 
when this maximum is reached, then from ^5j we obtain that the time interval 
that it takes for the function fi(t) to become zero is: 

tcollapse =1.1 (78) 

If we assume for instance that: 

/7 (t )=20xl0 3 kg/m 3 (79) 

we obtain: 

tcollapse = 465 S (80) 
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